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| ^ ■ Abstract. We introduce multi-sheeted versions of algebraic domains and 

quadrature domains, allowing them to be branched covering surfaces over the 
Riemann sphere. The two classes of domains turn out to be the same, and the 
main result states that the extended exponential transform of such a domain 
■ agrees, apart from some simple factors, with the extended elimination func- 

| tion for a generating pair of functions. In an example we discuss the algebraic 

curves associated to level curves of the Neumann oval, and determine which 
of these give rise to multi-sheeted algebraic domains. 
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1. Introduction 

Extending some previous work [23], [25], [26], [16], [17] on the rationality of the 
exponential transform we here go on to consider what we believe is the most 
general type of domains for which the exponential transform, in an extended 
sense, can be expected to have a "core" consisting of a rational function. Around 
this core there will then be "satellites" of some rather trivial factors, depending 
on the regimes at the locations of the independent variables. 

The "domains" we consider will actually be covering surfaces over the Riemann 
sphere. The terminology "quadrature Riemann surface" has already been in- 
troduced by M. Sakai [30J for the type of domains in question. An alternative 
name could be "multi-sheeted algebraic domain" , to extend a terminology used 
by A. Varchenko and P. Etingof [36]. In the present work we shall mostly use the 
latter terminology because we will not emphasize so much the quadrature prop- 
erties, but rather take as a starting point the way the domains are generated, 
namely by a pair of meromorphic function on a compact Riemann surface. 
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Given a bounded domain ficC, the traditional exponential transform [I], [24|, 
|25j . [26] , [T5] of is the function of two complex variables defined by 

1 f d( d( 
E n (z,w) = exp — / A 7 r 2,weC. 

For the unit disk D it is (see [15]): 

'1-^ MeC\D), 

(1) E(z,w) = 



1-4 (z6B,weC\: 

1-f (zGC\D,«)6 
G D). 



|z— u>| 



2 



1 



V 1— zw 

Using the Schwarz function [7j, [31] for d. 

S(z) = 
(pTJ) can be written 

(2) £M = (4^y {z) (^^y (w)s ( z ^ 

b{z)-w z-S{w) 
where £(z,w) is the rational function 

w) = = — 

zw 

and p = xd is the characteristic function of D. The expression ()2]) reveals the 
general structure of the exponential transform of any algebraic domain. 

The function S(z,w) is an instance of the elimination function, which can be 
defined by means of the meromorphic resultant TZ(f, g) of two meromorphic 
functions / and g. The resultant is defined as the multiplicative action of g 
on the divisor (/) of /, namely lZ(f,g) = g((f)), and the elimination function 
is £f >g (z,w) = lZ(f — z,g — w). In the case of the unit disk, or any quadra- 
ture domain, the relevant elimination function which enters into the exponential 
transform is the one with f(() = (, g(() = S((). This means that we in the 
present example get 

8{z, w) = n(C - *, 5(C) -*>) = (S(0 - tD)((C - z)) 
/^/>\ w, / n -i / \ \ S(z)—w zw — 1 

= (5(C) - W 1 • (Z - 1 ■ 00 = -Mr ~ = — , 

b{oo) — w zw 

as desired. 

The aim of the present paper is to generalize the formula ([2]) as far as possible. 
This will involve an extended exponential transform in four complex variables 
and an analogous extended elimination function in four variables, defined in 
terms of a conjugate pair of meromorphic functions on a fairly general compact 
symmetric Riemann surface. 
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The paper is organized as follows. Sections |2] and [3] contain general preliminary 
material. In Section 0] we introduce the concepts of multi-sheeted algebraic do- 
mains and quadrature Riemann surfaces and prove that they are equivalent. The 
main result is stated in Section [5] and proved in Section |6j Section [7J finally, is 
devoted to examples, namely the ellipse and Neumann's oval. 



2. The Cauchy and exponential transforms 

The Cauchy transform of a bounded density function p in C is 

c(z) = — MCRAdC 



2vri J C-z 

Typically the functions p which will appear in this paper will be like the char- 
acteristic function of a domain, or the corresponding integer valued counting 
function for a multi-sheeted domain. From C p , the density p can be recovered 
by 

(3) M = 

to be interpreted in the sense of distributions. 

If one writes the definition of the Cauchy transform as 

one realizes that it suffers from a certain lack of symmetry. A more balanced 
object would be the "double Cauchy transform" , 

(4) C p (z,w) = ± [p(0-^A 



27ri J C — z ( — w 
In fact, this double transform is much richer than the original transform, and 
after exponentiation it gives the by now quite well studied [4], [21], [2H], [2S], [15] 
(etc) exponential transform: 

E p (z, w) = exp C p (z, w). 
The original Cauchy transform can be recovered as 

C p (z) = res C p (z,w) = — lim wC p (z,w), 

at least if p vanishes in a neighborhood of infinity. One disadvantage with the 
double Cauchy transform is that the formula ([3]) turns into the more complicated 

dC p (z,w) = p(z) 
dz z — w 

On the other hand we have the somewhat nicer looking 

d 2 C p (z,w) 

-^— = -n P (z)5(z- W ), 
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where 5 denotes the Dirac distribution. 

Now, even the double Cauchy transform is not entirely complete. It contains the 
Cauchy kernel ^fj, which is a meromorphic differential on the Riemann sphere 
with a pole at £ = z, but it has also a pole at ( = oo. It is natural to make 
the latter pole visible and movable. That would have the additional advantage 
that one can avoid the two Cauchy kernels which appear in the definitions of 
the double Cauchy transform and the exponential transform to have coinciding 
poles (namely at infinity) . Thus we arrive naturally at the extended Cauchy and 
exponential transforms: 



(6) Cp (z,w,a } b) = ^- J p(C)( 



d( dc , A , d( d( 



C-z C-a (-w C~b' 
E p (z,w)E p (a,b) 



(7) E p (z } w;a,b) =expC p (z,w;a,b) 

Ep(z,b)Ep(a,w) 

If the points z, w, a, b are taken to be all distinct, then both transforms are well 
defined and finite for any bounded density function p on the Riemann sphere. 
For example, with p = 1, E p (z,w; a,b) turns out to be the modulus squared of 
the cross-ratio. See [16 1 for further details. 



3. The resultant and the elimination function 

Here we shall briefly review the definitions of the meromorphic resultant and 
the elimination function, as introduced in [16], referring to that paper for any 
details. If / is a meromorphic function on any compact Riemann surface M we 
denote by (/) its divisor of zeros and poles, symbolically (/) = / -1 (0) — / _1 (oo). 
If D is any divisor and g is a meromorphic function we denote by g(D) the 
multiplicative action of 5 on D. For example, if D = 1 • (a) + 1 • (b) — 2 • (c), 
a, b, c G M, then g(D) = ■ Now the meromorphic resultant between / and 

g is, by definition, 

n(f,g) = g((f)) 

whenever this makes sense. 
The elimination function is 

£fA z i w ) =K(f - z,g -w), 

where z, w € C are parameters. It is always a rational function in z and w, more 
precisely of the form 

where Q, P and R are polynomials, and it embodies the necessary (since M is 
compact) polynomial relationship between / and g: 

£ f MXO,g(0) = o (CeM). 
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We also have the extended elimination function, defined by 

./ - z g-w. 



£ f , g (z } w;a } b) = TZ( 



f -a g-b 



To relate the elimination function to the exponential transform one needs integral 
formulas for the elimination function. If / is meromorphic on M with divisor (/), 
let u f be a 1-chain such that dcrj = (/) and such that log / has a single-valued 
branch, which we denote Log/, in M\ supper/. Then Log / can be viewed as a 
distribution on M, and its exterior differential in the sense of distributions (or 
currents) is 



(9) dLogf = %-2mdH af . 

J 



Here dH af is the 1-form current supported by cr/ and defined locally, away from 
d(jf , as the differential (in the sense of currents) of that function H a , which is 
+1 on the right-hand side of cr/, zero on the left-hand side. Globally dH a is not 
exact (despite the notation), not even closed. To be precise, 

1 df 

d(dH af ) = — d(— ) = 5( f )dx A dy, 

Z7T1 J 

where 5(f) denotes the finite distribution of point masses (or charges) correspond- 
ing to (/). We shall also need the fact that dH a{ has the period reproducing 
property 



(10) / 6H aj A t 



M 



holding for any smooth 1-form r. 

Now we have (essentially Theorem 2 in [TB] ) 



v 11 ) £f, g (z,w,a,b) = exp[-^r f (-^ )ArfLog- — 



9 ~ v h 
b l 



The integrand is a 2-form current with support on the 1-chains <J g - w and <y g -b 
(because away from these curves the integrand contains d( A d(), so the integral 
is rather a line integral than an area integral. In fact, the above can also be 
written 

c ( 7 \ r f i d f df f df df 

E,„ {m b) = eMJ <J— Z - —J - j J— z - — a )l, 



which perhaps clarifies the connection to the definition of the elimination func- 
tion. 
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4. Multi-sheeted algebraic domains 

The boundary of a quadrature domain (algebraic domain) is an algebraic curve, 
but by no means every algebraic curve arises in this way. However, the gap 
between the two classes of objects can be reduced considerably by extending 
the notion of a quadrature domain, allowing it to have several sheets and to 
be branched over the Riemann sphere. This will take essentially one half of 
all algebraic curves into the framework of quadrature domains and exponential 
transforms. One step in this direction was taken in Sakai (30], where a notion 
of quadrature Riemann surface was introduced in a special case. Below we shall 
take some further steps. 

Let M be any symmetric compact (closed) Riemann surface. Slightly more 
generally, we shall allow M to be disconnected, namely to be a finite disjoint 
union of Riemann surfaces. The symmetry means that M is provided with an 
anticonformal involution J : M — > M, J o J = identity. If M is disconnected 
then J is allowed to map one component of M onto another. Let T denote the 
set of fixed points of J. Simple examples of symmetric Riemann surfaces are 
M = P (the Riemann sphere) with the involution being either Ji(C) = 1/C or 
-^(C) — ~~ 1/C 111 th e fi rs t case r = {( : |C| = 1} and M\T has two components, 
in the second case V is empty and hence M \ F has only one component. One 
can also think of identifying the points ( and J((). The identification spaces 
become, in the first case (J = J\) the unit disk together with its boundary, and 
in the second case (J = J2) the projective plane, thus a nonorientable surface. 

In general, the orbit space N = Mj J, obtained by identifying £ and J(() for any 
C G M, is a Klein surface, possibly with boundary. A Klein surface [2] is defined 
in the same way as a Riemann surface except that it is allowed to be nonorientable 
and that both holomorphic and antiholomorphic transition functions between 
coordinates are allowed. The possible boundary points of iV are those coming 
from T under the identification. From the Klein surface N, M can be recovered 
by a natural doubling procedure, described in Section 2.2 of [21] and in [2], 
for example. The latter reference actually describes several types of doubles 
(the complex double, the orienting double and the Schottky double), but the 
description in [31] will be enough for our purposes. In case iV is orientable and 
has a boundary the doubling procedure gives what is usually called the Schottky 
double, named after the inventor of the idea, F. Schottky [32J. The idea was 
later extended to more general surfaces by F. Klein [S]. 

Continuing the discussion of M, J and N, if M is connected but M \ V discon- 
nected, then M \ T has exactly two components, say M + and M_, and J maps 
each of them onto the other. It then follows that (M \ Y)/J can be identified 
with M + (or M_), and in particular that iV is orientable, hence is (after choice of 
orientation) an ordinary Riemann surface with boundary. If M \ V is connected 
then iV necessarily is nonorientable. 
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It is relevant to allow M to have several components. For example, in Section [7] 
we will encounter the double of the Riemann sphere, which simply is two Riemann 
spheres with the opposite conformal structure and with J mapping one onto the 
other. 

Now to the definition of "multi-sheeted algebraic domain". There are two in- 
gredients. The first is a compact symmetric Riemann surface (M, J) such that 
M is connected and such that M \ T has two components, one of which, call it 
M + , is to be selected. We could equally well have started with M + , to be any 
Riemann surface with boundary (bordered Riemann surface), and then let M be 
the double of M + . The second ingredient is a nonconstant meromorphic function 
/ on M. 

Definition 1. A multi-sheeted algebraic domainis (represented by) a pair (M + , /), 
where M + is a bordered Riemann surface and / is a nonconstant meromorphic 
function on the double M of M + . Two pairs, (M+, /) and (M+, /) are considered 
the same if there is a biholomorphic mapping : M — > M such that <poJ= J o 
and / = f o 0. 

The equivalence simply means that it is the image f(M + ), with appropriate 
multplicities, which counts. Note that also (M_,f) is a multi-sheeted algebraic 
domain, if (M+, /) is. Trivial examples of a multi-sheeted algebraic domain are 
obtained by taking M = P and J(Q = Then with / any nonconstant 

rational function (D, /) will be a multi-sheeted algebraic domain, as well as 
(P \ D, /). Some further examples will be discussed in Section [3 

Along with /, meromorphic on M, the symmetry J provides one more meromor- 
phic function on M, namely 

r = W^j)- 

With Q(z,w) the polynomial in (jS) for g = /*, the map M 9 (n> (/(C), /*(C)) 
parametrizes the curve Q(z,w) = (or, better, its projective counterpart). This 
parametrization is one-to-one (except for finitely many points) if and only if / 
and /* generate the field of meromorphic functions on M (form a primitive pair 
in the terminology of [9]). This will generically be the case if / is chosen "at 
random", but there are certainly many counterexamples. For example, / may 
be already symmetric in itself, i.e., / = /*, and then / and /* is a primitive pair 
only if M is the Riemann sphere and / is a Mobius transformation. If / and /* 
are not a primitive pair, then the polynomial Q(z,w) in (jSJ) is reducible. 

In case M + is planar (i.e., is topologically equivalent to a planar domain) and / 
is univalent on M + without poles on M + U T, then Q = f(M + ) is an ordinary 
algebraic domain, in other words a classical quadrature domain. In this case / 
and /* do form a primitive pair (see [13]). The reference to "quadrature" can be 
explained in terms of a residue calculation. In the present generality it is natural 
to use the spherical metric on the Riemann sphere in place of the customary 
Euclidean metric. This will remove some integrability problems, and the point of 
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infinity can be treated on the same footing as other points. Quadrature domains 
for the spherical measure have been previously discussed, at least in fluid dynamic 
contexts, for example in [36] (Hele-Shaw flow), [5] (vortex patches). 

Let h be a function holomorphic in a neighborhood of M + U T and assume, for 
simplicity, that / has no poles on dM + (such poles will actually cause no problems 
anyway). At least away from poles of / we have, by exterior differentiation, 



d 



hfdf \ hdfAdf 



l + ffj (l + ff) 2 
Therefore, if / has poles of orders rij at aj G M + then 

1 f hdfAdf , 1 f hdfAdf 
lim 



/M+n{|/|<l/e} (1 + ff) 
1 f hfdf 1 / hfdf 

— hm 



2vri J dM+ 1 + // e-o 2vri J m=1/e 1 + // 
1 f hf*df y- ? / \ \- hf*df 



Here it turns out that the terms Yl n jh( a j) cance l with corresponding terms 
(with negative sign) in the residue contribution unless /* happens to have zeros 
at the points aj. In any case, with or without cancellations, the right member 
above is of the form 

m «fc — 1 
k=l j=0 

i.e., equals the action on h by a distribution with support in finitely many points. 
The quadrature nodes bk are the poles of ■ 

If / is univalent the above identity becomes an ordinary quadrature identity 
(although for the spherical measure) of the form 



(12) l/ sW 

2tti J n (l + \z\ 2 y 



holding for all integrable analytic functions g in Q = f(M + ). Here the right 
member is given by L(g) = L(g o /), which still is the action of a distribution 
with finite support. 

When / is not univalent one should still think of the quadrature identity in the 
same way as in (fT2|) . just with the difference that Q is a region with several 
sheets over the Riemann sphere. The test functions should be allowed to take 
different values at points lying above one and the same point, but on different 
sheets. Thus those of the form g o f (i.e., those which would become g(z) in a 
formulation like ( TT2]) ) are too special. This is most easily expressed by pulling 
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everything back to M + , in which case we simply have the originally obtained 
identity 

(13) 27ri / M+ (TT17FF = 

for h holomorphic in a neighborhood of M + U T, and by approximation for func- 
tions h holomorphic and integrable (with respect to (i^/jf^ ) i n M + . An 
equivalent formulation is that there exists a positive divisor D in M + such that 



(14) ^1 h^y = , 



L [ h d ? Ad f 

2niJ M+ (l + |/| s 

holds for every holomorphic and integrable function h in M + with (h) > D. 

Definition 2. A quadrature Riemann surface (for the spherical metric) is a 
pair (M+, /), where M + is a bordered Riemann surface and / is a nonconstant 
meromorphic function on M + such that J M+ (^jjf^ < oo and such that (|T3|) 
(or (THj) ) holds for some L (respectively -D) and the indicated classes of functions 
h. 

The notion of equivalence between pairs is the same as in Definition [TJ If / is 
meromorphic on M then we have 

If df Adf If df Adf n f dz A dz 

— / n , .x.oxo = 7T-7 / , ...mo = n < OO, 



27ri 7 M+ (1 + |/| 2 ) 2 " 2vri J M (1 + |/| 2 ) 2 2vri J F (1 + |^| 2 ) 2 
where n is the order of /. Thus we have proved one direction of the following. 

Proposition 1. A pair (M+, /) is a multi-sheeted algebraic domain if and only 
if it is a quadrature Riemann surface. 

Proof. It remains to prove that a quadrature Riemann surface is a multi-sheeted 
algebraic domain. In many special cases this has already been done (for the 
Euclidean metric), see for example [I], [13], [30]. We shall discuss here the general 
case under the simplifying assumption that / is meromorphic in a neighborhood 

ofM+ur. 

It will be convenient to introduce some further notation. With A G M any subset 
and D any divisor on M we denote by Od(A) the set of functions h meromorphic 
in a neighborhood of A and with (h) > D. Similarly, Ojf(A) denotes the set 
of 1-forms u, meromorphic in a neighborhood of A and satisfying (u) > D. We 
shall also use standard notations for cohomology groups. 

Now, returning to the previous residue calculation one realizes that what needs 
to be proven is that if 

(15) ±f A ^L = 

2xii 81 . 1 + // 
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holds for all h G Od(M + UT), for some sufficiently large divisor D, then / extends 
to a meromorphic function on M. On T we have 

fdf f*df 



1 + ff 1 + //*' 

where the right member is holomorphic in a neighborhood of T and can be 
viewed as representing an element in the cohomology group H\M,O]l D ). When 
D is strictly positive this group is trivial since, by Serre duality [33], [TO] . 
H\M,d%) H°(M,O d )* = O d (M)* = (here star * denotes dual space), 
hence there exist u)± G Ojjj(M± U Y) such that 

f*df 



in a neighborhood of T. 
Clearly 



1 + //* 



1 

2tti 



uo + — w_ 



htu + = 

'<9M+ 

for h G (9i)(M + ur) since the integrand is holomorphic in M+, hence (fT5l) reduces 
to the statement that 

(16) -L f hu„ = 

2m JdM+ 

for h G (9d(^+ U r). At this point we may use a general duality theorem (also 
related to Serre duality) going back to the work of J. Silva [35] and G. Kothe [18], 
extended by A. Grothendieck [12j and, in the form we need it, C. Auderset [3]. 
It states that the bilinear form in h and Co>_ defined by the left member of (I16p 
induces (with the path of integration moved slightly into M_) a non-degenerate 
pairing 

O d (M + U T)/O d {M) x Oh D {M.)/Oh° D {M) -> C, 
which exhibits each of the quotient spaces as the dual space of the other. 
In view of this ffTB"]) implies that w_ G O 1 J_ D (M), in particular that 

(17) = oj + - c_ G o^(M + u r). 

Since / is meromorphic in a neighborhood of M + U T, (JTTj) implies that also /* 
is meromorphic there, hence that / actually is meromorphic on all M, as was to 
be proved. 

The usage of the general duality theorems above can be replaced by more direct 
arguments, like applying ( TT3"]) to suitable Cauchy kernels. Specifically we may 
choose, with ( G M_, 

h(z) = $(z,C;^o,CoR, 
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where the right member is a kernel which in the case of the Riemann sphere is 
the usual Cauchy kernel 

*(*,C;*b,Co)dC = 7^-7^-, 
C - z Q - zo 

and which has counterparts with good enough properties on all compact Riemann 
surfaces (see [27] ). The point Co is needed in higher genus. Actually the duality 
theorem discussed above can be proved using this kernel. 



5. Statement of the main result 

Let M, M + , T = <9M + , / be as in Section HI more precisely such that (M + , f) 
is a multi-sheeted algebraic domain. To account for the multiplicities of f(M + ) 
as a covering of the Riemann sphere we introduce the integer-valued counting 
function, or mapping degree, 

p(z) = card {(GM + : /(C) = z}, 

pointwise well-defined for z G P\/(r). It is understood that points C are counted 
with the appropriate multiplicities. Set also 



(18) S = /oJo/-i = /'o/- 1 . 

This is a multi-valued algebraic function in the complex plane which contains all 
local Schwarz functions of f(T), because for z G f(T) one of the values of S(z) 
is z. 

We will have to make expressions like (S(z) — w)) p ^ well-defined, i.e., single- 
valued, despite S(z) itself being multi-valued. The natural definition is the fol- 
lowing: 

(s(z)-wy^ = (r-w)((f-z)\M + ). 

Here (/ — z)\m + denotes the restriction of the divisor (/ — z) to M + and the right 
member then is the multiplicative action of /* — w on (/ — z)\m + - To spell it 
out, let 

f- 1 (z)riM + = {( 1 ,...,( P{z) }, 

with repetitions according to multiplicities. Then 

(19) (S(z) - w)M = (f (Ci) - w) (/*(Cp«) - w), 

which is a natural definition in view of ffl8l) . Clearly (S(z) —w) p ^ is an analytic 
function of z in regions where p(z) is constant. 

Now, for the main result, we have two functions which we want to relate to each 
other: one is the weighted exponential transform 



E p (z, w; a, b) = exp[^r jT p(C) ^ 



^C cK_\ /_d£ dty 

C-z (-a) \Q-w C-K 
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^iJM+Kf-z f-aj \f-w f-b, 
which can be viewed as a kind of potential of p, and the other is the elimination 
function, which is defined by algebraic means and always is a rational function, 
namely of the form 

e ( - ts vff~ z f*-w. Q(z,w)Q(a,b) 
E U . u, a, 6) = K(— a , 77 - T ) = Qizrb)Q(aM - 

The latter expression comes form (jSj) together with the observation that the one 
variable polynomials cancel in the four variable case. The nature of E p (z, w, a, b) 
depends on the locations of the points z, w, a, b, more precisely on the the values 
of p at these points. The main result is the following. 

Theorem 2. Let (M + ,f) be a multi-sheeted algebraic domain. Then, in the 
above notations, 

E p (z, w; a, b) = £fj*(z, w; a, b)- 
z — w \ w — z \ I a — o x ' — " x 



S(z)-w) ys(w)-zj \S(a)-bJ \S(b) - a 
S(z)-b\ p(z) (s{^)-aY {W) fS(a)~w\ p{a) f ~S(bj - ^ ^ 



z — b ) \ w — a J \ a — w 

6. Proof of the main result 

Extending (j^D to four variables gives 

dC p (z,w;a,b) _ p(z) p{z) 



dz z — w z — 6 ' 

which tells that the function 

z — b 

C p (z, w; a, b) + p(z) log 

z — w 

is analytic in z in regions where p(z) is constant, namely in each component of 
P \ /(r). Hence so is the exponential of it, namely 

E p (z,w;a,b)-(^—ty^. 

z — w 

Augmenting this we have that also the function 

„T.\ P( Z ) 



F(z) = E p (z, w; a, b) 



z — b S(z) — w 



z — w S(z) — b 

is analytic in P \ f(T), away from poles caused by the presence of S(z). Now we 
claim that it is even better than that: F(z) is meromorphic everywhere, hence is 
a rational function. To prove this it is enough to prove that F(z) is continuous 
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across f(T). It is well-known that E p (z, w; a, b) is continuous in z (see for example 
[T5]). so we only have to bother about the other factor. But it is easy to realize 
that also this is continuous: spelling out as in ( Fl9|) and assuming for example 
that p(z) increases by one unit as f(T) is crossed at a certain place we find that 
the factor 

'S(z)-w\ p{z) 



w 



changes from 



(/*(&) -a) (r(Cpw)-w 



to 

(/*(Ci) - g) (/%(*)) - w)(r(( p(z)+1 ) - w) 

(/(CO - ™) (/(Cp (2) ) - <d)(/(Cpm +1 ) - ' 

which clearly is a continuous change since the new point Cp(*)+i starts up on T. 
Similarly for the factor 



Repeating the above argument for w, a, b it follows that the function 



, z — b S(z) — w\ I w — a S(w) 
E p (z,w;a,b) 1 



z — w 



p(w) 



S(z)-bJ \w-z S(w) 



a-w S(a)-b\ pia) fb-z S(b)-^ P(b) 



— b S(a) — w J \ b — a S(b) — 

is rational in the variables z w, a, b. Thus, since also £fj*(z,w;a,b) is rational 
in these variables, it is enough to prove that the formula in the statement of the 
theorem holds just locally, somewhere. We may then choose z, w, a, b close to 
each other, so that in particular p(z) = p{w) = p(a) = p(b). In addition we 
may assume that this value is the smallest value of p occurring on P. The case 
that it is zero can be treated exactly as in the proof of Theorem 6 in [16] , which 
concerns the special case a = b = oo. 

So let us for example assume that p(z) = p(w) = p(a) = p{b) = 1 (the general 
case will be similar). Thus / attains the values z, w, a, b exactly once in M + , 
and these four points on the Riemann sphere are close to each other. Let 7 be an 
arc in P from b to w (e.g., the geodesic arc). Then the function z 1— > log has 

a single- valued branch, call it Log in P\ 7, hence Log is single- valued in 

Z O J z 

M \ J -1 (7). We may consider / _1 (7) as a 1-chain, and as such it has the same 
role for Log^£^ as 07 has for Log / in so that 

/ - W d f d f o • , TT 
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If / has degree n, then / ^7) consists of n small arcs, one of which is located 
on M + . Let o = / _1 (7) H M + be that arc and let a = J (a) be the reflected arc 



in M_. In the sequel we shall use / 1 in the restricted sense / 1 = (/|m 4 
Thus 



da = r i ( W )-r i (b) 

da = FHw) - fHb) 



f _ 5 ' IM +' 
7* -6 



M_ 



Using dTTJ), (J9J) and (JIDj) we now get 

1 f t df df w JT „ f-^ 



£ fJ *(z,w;a,b) = exp[-^ / (-^ )AdLog 

27T1 Jm J - % J ~ a 



2vri J M+ /-z /-a /* - 6 J* f-z f-a 

Here we start by rewriting the last factor according to 

df df f - z f^,) 



exp / ( 

J a 



zf-z f-a' l f-a s f-Hb) 



/*(/-») ~ g . Hf- l {b)) - a = SW) - g S( W )) - a 
/*(/-») - a ' /*(/- 1 (6)) - * SH) - a ' - z' 



Note that functions like S(w) — z = (S(w) — z) p ( w ^ are single-valued in the 
present case (cf. (EES}). 

Next, the first factor can be integrated partially, to become 

1 f , df df x A JT /*-<z> 

exp 



2 ^Jm+ f-z f-a 

df df 



2m JdM+ 



f-z f-a' 



1 f , ( df df 

' exp — / d (- -z Log 

2tti J m+ f-z f-a 

6XP[ 2ri] dM }f-z f-a' 



Mi 



6 /* 


- 6 


r 


— w 


/* 


-b 


r- 








/* 


— w 


,g /* 


-b 


/*- 


- w. 



1 f df df f*-w, S(z)-w S(a)-b 



2niJ dM f-z f-a f*-b S(z)-b S(a)-w 
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We have to rework all expressions which are not yet in the form appearing in 
the statement of the theorem. So we next turn our attention to the first factor 
in the last obtained expression. This can be rewritten as 

r 1 f ( d f df f*-w 1 



2vn J dM+ f - z f-a f* -b 

r 1 

exp[ 



2m 



dM 4 



«q>[~ / d 



M 



+ 



df 


df 


)Log 


./- 


w 


f-z 


f-a 


' f- 


b 


df 


df 


)Log 


f- 


w 


■f-z 


f-a 


' /- 


b 



r 1 f , df df f~w 1 
ex P[7T^ b 7 )dLog- 



27riJ M+ f-z f-a f-b 
= exp[- f (Sf-i (z ) -S f -i(a)) dxdy Log ^ W - } 

J M+ J — V 



2niJ M+ \f-z f-a) \f-w f-b 

f ( df df — 

•exp / (- )dH a 

Jm+ f-z f-a 

z-b a-w 1 f ( df df \ f df df 
■ exp[— / A 



z-w a-b 2m J M+ \f - z f-a) \f - w f-b 

• ex P - / 7 7 

J a f-z f-a 

z — b d — w w — a b — a 

rjp{z, w; a, b) 



z — w a — b w — z 



Now putting all the pieces together we obtain the formula in the statement of 
the theorem. 



7. Examples 

7.1. The ellipse. Let D denote the domain inside the ellipse 

— + yl = i 

a 2 b 2 

with semiaxes a > b > and foci ±c = ±V a 2 — b 2 . The exterior domain ¥\D 
is known to be a null quadrature domain for the Euclidean metric. For 
the spherical metric it is a two point quadrature domain. Indeed, the Schwarz 
function for the ellipse is 

S(z) = a ^z± 2 -^Vz^ 2 , 
c 2 c 
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and for h holomorphic in P \ D and smooth up to the boundary we have 

1 f dz A dz v-^ S(z)dz 

tt^ I h{z)— . = > res h(z)- tttt- 

2wi J A d KJ (l + \z\ 2 ) 2 ^ z &\d y J l + zS(z) 

The residues come from the zeros of 1 + zS(z) in the exterior of the ellipse, 
and straight-forward computations show that there are exactly two such zeros, 
located on the imaginary axis and more precisely given by 



1 



z = ±z = ±— V a 2 + b 2 + 2a 2 b 2 + 2abVl + a 2 + b 2 + a 2 b 2 . 
ic v 

Thus we have a quadrature identity of the form 

(20) JL J^ h(z) JEA^_ = Co(ft(2o) + h( . :o)h 

c being the residue of f^f^ at z = ±z . 

The exterior of the ellipse is the conformal image of the unit disk under the 
Joukowski map 

c 2 C 2 + (a + bf 



/(C) 



2(a + 6)C 



Thus (D, /), or simply (P \ D, z) with z denoting the identity function, is a 
(single-sheeted) algebraic domain. The same function / maps the exterior of the 
unit disk onto a multi-sheeted algebraic domain, i.e., (P\D, /) is (or represents) 
such a domain. It covers D twice and P \ D once, in other words the counting 
function is 



2 for z E D, 
1 for z G P \ D. 



Since there are several sheets the associated quadrature identity is best expressed 
in a form pulled-back to P \ D, i.e., on the form (Tl3|) . A slightly weaker form is 



^- [ p(z)h(z) Z A I = ci(/t(zi) + h(-z 1 )) 1 

Z7T1 Jp [L -+- \Z\ ) 



where 



ic 



±zx = ±-ya 2 + b 2 + 2a 2 b 2 - 2abVl + a 2 + b 2 + a 2 b 2 



and c\ is the residue of x+zS(z) a ^ z = ^ z i- This form is weaker because it only 
uses test functions h(z) that take the same values on the two sheets over D. 
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7.2. Neumann's oval: classification of level curves. By inversion in the 
unit circle and a rotation by 90 degrees (for convenience) the ellipse transforms 
into a curve known as Neumann's oval [22], [23], [31], [19] with equation 

(21) a 2 b 2 (x 2 + y 2 ) 2 -a 2 x 2 -b 2 y 2 = 0. 

The exterior of the ellipse transforms into a bounded domain Q, and since in- 
versions and rotations are rigid transformations with respect to the spherical 
measure, also Q will be a two point quadrature domain for the spherical mea- 
sure. The formula is immediately obtained by inversion and rotation of f l20|) . The 
domain Q also satisfies a quadrature identity for the Euclidean measure (indeed, 
both types of quadrature identities are equivalent to S(z) being meromorphic in 
Q). The latter quadrature identity is somewhat simpler, namely 

It turns out to be a quite rewarding task to investigate the algebraic curves 
corresponding to the level curves of the left member in (I2ip , and in particular to 
determine which of them correspond to multi-sheeted algebraic domains. Most 
types of phenomena which could possibly show up really do show up among these 
curves. This task is what we are going to undertake for the remainder of this 
section. 

In order to simplify a little we first scale so that the quadrature nodes above 
become ±1. This means that 2ab = c. Then we need only one parameter (in 
place of the two, a and b), which we take to be 

vV + b 2 
r = > 1. 

c 

The quadrature identity now becomes 

(22) / hdxdy = 7rr 2 (h(-l) + h{l)), 

Jn 

holding for all integrable analytic functions h in Q. Set 

-O/ \ 2 2 2 2 o 2 

Q{z,w) = z w —z — u> — 2r zw, 

which is the polarized version ((z,z) polarizes into (z,w)) of the left member 
in ( 12T|) . There are exactly two open sets for which the quadrature identity (|22|) 
holds, namely 

Q = {z e C : Q(z,z) < 0} 

and 

[n] = {z e C : Q(z,z) < 0} U {0}. 

The latter is just the completion of the former with respect to one missing point. 
The domain Q (or [Q]) can be viewed as two disks glued together, or "smashed", 
or "added", and has been studied by many authors, for example [6], [19], [20] . 
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To analyse the level curves of Q(z, z) we set, for any «6l, 
(23) Qa{ z i w ) — <2 2 w 2 — z 2 — w 2 — 2r 2 zw — a. 

Let 



Q a (t, z, w) = z 2 w 2 - z 2 t 2 - w 2 t 2 - 1r 2 zwt 2 - at 4 



be the corresponding homogenous polynomial. We shall keep r > 1 fixed and 
just vary a. The real locus of Q a in C is 



It represents the intersection with {w = z} ("the real") of the complex locus in 



loc c Q a = {(z, w) eC 2 : Q(z, w) = a}, 
which has a natural completion in the projective space Pa(C) as 

loc Q a = loc P2{c) Q a = {(t : z : w) E P 2 (C) : Q a (t, z, w) = 0}. 

Here C 2 is embedded in F 2 (C) so that (z,w) corresponds to (1 : z : w). By 
real points in C 2 we mean points (z, w) satisfying w = z, and these are the 
fixed points of the involution J : (z,w) h-> (w,z). In projective coordinates the 
involution is J : (t : z : w) i— y (t : w : z). In addition to the anticonformal 
involution J, the curve loccQ a has the conformal symmetries (z,w) h-> (w,z) 
and (z,w) h-> (—z,—w). 

Let (M a , J a ) be the compact symmetric Riemann surface corresponding to (locQ a , J). 
As point sets they are identical except for a few singular points on \ocQ a , which 
are resolved on M a . The real locus of Q a corresponds to the fixed point set r a 
of J a . 

A first observation is that the function z i— > Q(z, z) has five stationary points in 
the complex plane. There are two global minima, on the level a = — (r 2 + l) 2 , 
there are two saddle points on the level a = — (r 2 — l) 2 and there is one local 
maximum (at the origin), on the level a = 0. These three values of a will 
correspond to changes of regime for the algebraic curve Q a (z,w) = 0. 

Solving the equation Q a (z,w) = for to as a function of z gives the Schwarz 
functions for the curves in the real locus. The result is 



We see that S a (z) in general has four branch points. The levels a at which Q(z, z) 
has stationary points are exactly those values of a for which some or all of these 
branch points resolve: for a = — (r 2 ± l) 2 the square root resolve completely into 
second order polynomials, and for a = one pair of branch points resolves. 

Let p = p(a) denote the genus of M a , or equivalently of locQ a . The degree of 
Q a is four, hence the genus formula in algebraic geometry [TT] , [2T] tells that 



loc R Q a = {z E C : Q(z, z) = a}. 



C 2 



(24) 




p + s 



(4-1). (4-2) 



3 



2 
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where s > is a certain number related to the singular points. An analysis, 
carried out in detail in [TJ], shows that \ocQ a passes through the points (0 : 
1 : 0) and (0:0:1) and that it at each of these points has two simple cusps 
of multiplicity one with distinct tangent directions. In particular, the points 
(0:1:0) and (0:0:1) are singular, and it turns out that each of them gives the 
contribution +1 to s. Except for the above two points of infinity, locQ a stays in 
C 2 . Thus, what remains of the genus formula is 

(25) p + e = l, 

where e denotes the contribution to s which comes from finite singular points. 
By ([2"5]) . e < 1, so there is at most one finite singular point, and this must be 
visible in the real because nonreal singular points necessarily come in pairs. 

The above analysis preassumes that the curve locQ a , or polynomial Q a , is irre- 
ducible. This is the case for most values of a, but there are two exceptions: 

(i) For a = — (r 2 — l) 2 , Q a factors as 

Q a (z, w) = ((z + l)(w + 1) - r 2 ){{z - l)(w - 1) - r 2 ). 

Each factor defines its own algebraic curve and Riemann surface. These have 
genus zero and are moreover symmetric: the zero locus of each factor is preserved 
under the involution (z,w) (w,z). The real locus \oc^Q a is the union of two 
intersecting circles, those of radius r and centers ±1. The intersection points, 
z = ±\/r 2 — 1 are saddle points for the function z i— >■ Q(z, z). 
{ii) For a = — (r 2 + l) 2 , Q a factors as 

Q a (z, w) = {{z + l)(w - 1) - r 2 )((z - l)(w + 1) - r 2 ), 

where again each factor defines its own algebraic curve and Riemann surface of 
genus zero. However, in the present case they are not symmetric, instead the 
involution maps each of these Riemann surfaces onto the other. The real locus 
consists only of the two points ±r. This can easily be understood by observing 
that the value a = — (r 2 + l) 2 is the global infimum of Q(z, z). 

In both of the reducible cases Bezout's theorem says that there should be four 
points of intersection between the two curves (since these have degree two). These 
intersection points are the two points of infinity (0:1:0) and (0:0:1) plus, in 
the first case the intersection points of the two circles (in the real), and in the 
second case the two points ±r. 

Besides the above two special values of a, also the quadrature value a = is 
exceptional. This is a local maximum value for Q(z, z). The local maximum is 
attained at z — 0, and the corresponding point (0, 0) on the algebraic curve is a 
singular point (since both partial derivatives of Q vanish there). Thus e = 1 in 
the genus formula f )25|) . hence p = 0. 

We now embark the full classification. Pictures for the case r = a/2 with a = 
3,-0.5,-1,-1.5 are shown in figures [1] and [21 where the shaded areas are the 
sets where Q a (z, z) < 0. 
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For a > 0, T a has one component and there are no singular points visible 
in the real. Therefore, by (|25|) and the remark following it M a has genus 
one. This means that the symmetry line r a is not able to separate M a 
into two halves (see more precisely discussions in Section 2.2 in [31]). Thus 
M a \ T a has only one component, and it will not generate any algebraic 
domain (even multi-sheeted), despite the nice picture in the real, with a 
smooth algebraic curve bounding a simply connected region (figure [U left). 
M a can be viewed as the double of a Mobius band. 

At a = the genus of M a collapses to zero, T a has still only one compo- 
nent even though IocrQc* has gotten an additional point, (0,0), which is 
a singular point of locQ a (resolved on M a ). It follows that M a \ T a has 
two components, one of which, say M + , is mapped conformally onto the 
quadrature domain [Q] by the analytic function / which corresponds to the 
projection (z, w) (-)■ z on locQ a . 

For — (r 2 — l) 2 < a < the genus of M a is again one, and the singular 
point (0, 0) in the previous case has now grown up to a curve, hence T a has 
two components. Also M a \ T a has two components, say M±, and M a can 
be viewed as the double of M + (or M_), which topologically is an annulus. 

The meromorphic function / : M a — > P which corresponds to (z, w) i— >■ z 
on locQ a is however no longer univalent (not even locally univalent) on 
what corresponds to M + in the previous case. Therefore this gives a now 
only multi-sheeted algebraic domain, with f(M + ) consisting of a main piece 
which contains the points ±1 and the origin (the dashed area in figure [TJ 
right) plus a smaller piece around the origin (the bounded undashed region). 
The latter piece thus is covered twice, and the two sheets are connected via 
two branch points of the Schwarz function. 

For a = — (r 2 — l) 2 the curve is reducible, hence M a is the union of two 
independent Riemann surfaces, both of genus zero and symmetric under 
J a . In the real locus 1ocrQ q we simply have two intersecting circles, those 
centered at ±1 and having radius r (figure [21 left). Explicitly: 

Q Q (^) = (|^-l| 2 -r 2 )(|^ + l| 2 -r 2 ). 

M a is the union of two Riemann spheres and can be viewed as the double 
of two disks. 

For — (r 2 + l) 2 < a < — (r 2 — l) 2 we are back to the case of genus one with 
T a and M a \ T a both having two components. However, the situation has 
changed in the sense that the involution goes the other way (like z > — z 
in place of z t— > z in a right-angled period parallelogram), it may be more 
natural in this case to think of M a as the double of a cylinder than as 
the double of an annulus (even though these two types of domains are 
topologically equivalent). The real locus consists of two closed curves, one 
enclosing two branch points close to z = 1, the other enclosing two branch 
points close to z = — 1, and with / : M a — > P as before, / maps M + (say) 
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onto the dashed region to the right in figure [2] (right) covered twice and the 
(unbounded) undashed region covered once. 

• When a = — (r 2 + 1) 2 the two closed curves in the real locus of the previous 
case have shrunk to two points, the minimum points z = ±1 of Q(z,z). 
The curve is reducible and M a hence is the disjoint union of two Riemann 
surfaces (of genus zero), but these are not symmetric under J. Instead 
the involution J maps each of them onto the other. The algebraic curve 
IockQq consists of two pieces, which meet each other in two points of 
tangency. This is what is seen in the real locus. Thus M a can be viewed 
as the double of the Riemann sphere, and J a has no fixed points. 

• When a < — (r 2 + l) 2 finally, there is no real locus at all. The genus is 
one, but now J has no fixed points at all (r Q is empty). Therefore no 
algebraic domain (even multisheeted) can be associated to this case. M a 
can be viewed as the double of Klein's bottle (see again [2J, and also [3T] 
for doubles of nonorientable surfaces). 

We summarize the discussions as follows. 

Proposition 3. The surface M a (possibly disconnected) can be viewed as the 
double of the Klein surface N a = M a /J a , which in the different regimes ofa&M. 
is of the following topological type. 

• For a > 0: a Mobius band 

• For a — : a disk. 

• For — (r 2 — 1) < a < 0: an annulus. 

• For a = — (r 2 — l) 2 : two disjoint disks. 

• For — (r 2 + l) 2 < a < — (r 2 — l) 2 : a cylinder. 

• For a = — (r 2 + l) 2 : a sphere. 

• For a < — (r 2 + l) 2 : a Klein's bottle. 

The pair (M + , /) defines a single- sheeted algebraic domain for a = and a multi- 
sheeted algebraic domain for — (r 2 — l) 2 < a < and — (r 2 + l) 2 < a < — (r 2 — l) 2 . 
For a = — (r 2 — l) 2 it defines two single- sheeted algebraic domains (which intersect 
in the complex plane). 

The exponential transform of the above multi-sheeted algebraic domains is given 
by Theorem [H where the elimination function is 

S fJ . (z, w; a, b) = — — 
Q a {z,b)Q a {a,w) 

with Q a as in (23) and the Schwarz function is given by [Efy . 
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Figure 1. The level set T a for r = a/2 and for a = 3 and a = -0.5 
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Figure 2. The level set r Q for a = -1 and a = -1.5 
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